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Abstract 

We present a microscopic theory for both equilibrium and nonequilibrium 
transport properties of coupled double quantum dots (DQD). A general for- 
mula for current tunneling through the DQD is derived by the nonequilibrium 
Green's function method. Using a Hartree-Fock approach, effects of multi- 
level coupling and nonequilibrium electron distributions in resonant tunnel- 
ing are considered. We find that the peak in the resonant tunneling current 
through two symmetric dots will split only when the inter-dot coupling is 
stronger than dot-lead coupling. We predict that population inversion can be 

achieved in one dot in the nonequilibrium regime. 
73.20.Dx, 73.40.Gk, 73.50.Fq 
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Coupled quantum dot systems have received much attention recently |T]-§]. Resonant 
tunneling through zero dimensional (OD) states of coupled quantum dots has been studied 
experimentally by several groups in both equilibrium "C ksT) and nonequilibrium 

(\fiL— Hr.\ ksT) regimes (where fii and are chemical potentials of leads externally 
attached to the quantum dots from left and right, respectively). The Coulomb blockade 
theory (CBT) was used to study the equilibrium properties of electronic tunneling in a 
double quantum dot (DQD) system and more recently [§J§, has been extended to explain 
the tunneling current peak splitting observed in the experiment of Waugh et al |§ . 

Most previous theoretical studies on DQD systems @-§|J^Q] concentrated on the equilib- 
rium properties. Interesting properties in the nonequilibrium regimes @, where each dot can 
have different thermal nonequilibrium states driven by the two leads, have received much less 
attention. In this paper, we use a microscopic tunneling model to study the electronic trans- 
port through DQD systems based on the nonequilibrium Green's function (NGF) method. 
Using the NGF method, the resonant tunneling current can then be derived exactly in both 
equilibrium and nonequilibrium regimes for tunneling through an interacting quantum dot 
|TT| and non-interacting multi-quantum dots [ XT | . Here we give a general expression for 



the current for tunneling through interacting multi-level double quantum dot systems. The 
general current we derived can be used to obtain the well-known current equations in the 
equilibrium regime. In the nonequilibrium regime, it is usually difficult to satisfy detailed 
current balance between different dots and leads in various approximations ||12|| , except in 



the case that the contribution to the imaginary self-energy is only from the lead-dot cou- 
pling. In this case, the current equation is much simplified and can be written in a compact 
form. 

In this paper, we study resonant tunnelings through double quantum dots, each of them 
having more than 50 electrons. In these systems, only the quasi-particle states are relevant 
to the resonant tunneling. The spectrum of quasi-particle states near the Fermi surface will 
not change significantly with changing the electron number from N to N + 1 except for a 
slight energy level shift mainly due to the "charging effect". This has been confirmed for 



quantum dot systems with N > 30 by microscopic calculations |[20|| . In general, various scat- 
tering contributions to the imaginary self-energy (scattering-in scattering-out processes) can 
be discussed phenomenologically using our current equation below, which is useful for inter- 
pretation of experiments. We will use the general current equation derived here to discuss 
experimental results performed in the nonequilibrium regime 0; some of our explanations 
are different from previous ones ||. In particular, we will see that the study of nonequilib- 
rium electron distributions in the quantum dots are crucial to understand nonequilibrium 
experiments [[|. If the couplings of lead-dot and dot-dot are T and \t\, respectively, we pre- 
dict that in symmetric equilibrium experiments H a peak of the resonant tunneling current 
will split only when \t\ > T/2 with a splitting magnitude 2(|t| 2 — r 2 /4) 1//2 . 

Particularly, we predict that a population inversion can be achieved in one of the two 
quantum dots when the tunneling between the two dots is in resonance. This population 
inversion should be easily realized under conditions such as the low temperature experiments 
of Ref. 0,0]. At high temperatures, the population inversion critically depends on electron- 
phonon (e-p) or electron-electron (e-e) scattering induced relaxation times. This population 
inversion is similar to that in quantum cascade lasers (QCL) [[DJ except that the former 
co-exists with resonant tunneling and the tunneling is not photon-assisted. Thus population 
inversion can be achieved in the double quantum wells under similar conditions as for QCL. 

The structure of this paper is as follows: in section | we discuss the model and derive the 
general current equations; in section we discuss various approximations and the simplified 



current equations; in section III we discuss the non-equilibrium distribution of electrons in 



a double quantum dot system and the possibility of population inversion; in section [IV|, we 
discuss two recent experiments; and we conclude in section |V|. 

I. MODEL AND CURRENT EQUATIONS 



We model the coupled DQD system with two attached leads by the Hamiltonian [13 
H = H + H t + H h with 
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i=l,2 



(1) 



where the subscripts (L, i?) represent the left- and right-lead, and i sums over the two 
dots. In Hartree-Fock (HF) representation, ipl G{L , R} = (4 lQ , a\ 2Cl , ■ ■ ■ , aj^J, ^J e{1)2} = 
( c !,i> c l,2> ' ' ' 5 c !,nJ- Here {ki, ■ ■ ■ k^ a } and {1, 2, • • -rii} are labels of states in the lead and 
dot, N a and are the numbers of total states in the lead-a and dot-i, respectively. The 
E*'s are diagonal matrices whose elements are calculated in the HFA. H t is a tunneling 
Hamiltonian. V a ,i = [^f'M i s the (N a x m) hopping integral matrix; the matrix element 
V^i is the hopping integral from state k m in lead-a to state-/ in dot-z. t i2 = [t; 1 '^] is a 
(ni x n 2 ) matrix, and the matrix element t]^ is the hopping integral from state-/ in dot-1 to 
state-m in dot-2. Hj describes residual interactions such as electron-phonon (e-p) scattering, 
disorder scattering, and e-e correlations. 

Since the system under investigation is in non-thermal equilibrium in the nonlinear region 
(with finite bias voltage), The natural technique to use is that of non-equilibrium Green 
functions [pl6|-|i~8|j (NGF). The NGF in the DQD in matrix form is defined by: 



where Tc is the chronological time-ordering operator on the contour C of the closed time 
path, and r = t) with j3 = ± in the plus (minus) branch of the contour C. We choose 
the bases as HF states, so the evolution of ipi is governed by H t and Hj. Similarly, we define 
the GF for the two leads Gl and Gr, and the "off-diagonal" NGF G^i and G^, etc. 
The four combinations of the time branches (±, ±) give four linearly dependent component 
NGFs. Here we have used the notation that G 7i< 5 without superscript index means a (2 x 2) 
matrix of G^f with superscripts. After a rotation in the (2 x 2) space, we can obtain 
three independent components of the NGF, i.e., G[j a ' + ~ (retarded, advanced and statistical 
component carrying information of nonequilibrium distributions). Another frequently used 
(Keldysh) statistical component of NGF instead of G +_ is G K = 2G +_ + G r - G a . Each 




(2) 
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of the three components is itself a (r^ x rij) matrix, where n« is the number of total states 
in dot-i or lead-i. (Note that we have used the notation rij for the leads instead of N a , this 
is only for the convenience of description here). The matrix elements of these rii x rij NGF 
matrices are just the scalar NGF, e.g. [G[ 2 ]i,mit,t') = — iQ(t, t')({c n (t), cj^^t')}}. 

Since the leads are open systems, we can treat them as two equilibrium Fermi seas of 
non-interacting quasi-particles with different chemical potentials. The degrees-of-freedom of 
the leads can be integrated out and the self-energies (in matrix form) of the two quantum 
dots due to the tunneling between the leads and dots become S°(o;) = ct^Vl GlV-lcTx] 
S^tu) = ct x Vr^GrVrct x , where a x is the Pauli matrix operating only in the (2 x 2) Keldysh 
space g7|. Defining iTj = (E°) r - (£°) a , it is easy to see that: (E°) + -(w) = if L {w)Ti(u) 



and (Y>2) + ~(uj) = if R (u)T 2 (uj), with /l(u;) and /r(u;) the Fermi-Dirac distribution functions 
of the left and right leads, respectively. Physically, the diagonal matrix element [Fj]^ of Tj is 
the level broadening of state-/ in dot-?, due to the tunneling to all the states in the neighboring 
lead. 

The current equations can be derived formally using the Dyson equations. In our system, 
the Dyson equations (in Keldysh space) |HJ can be written as: 

Si Gy = 5{j + SjjGjj + a x t\2Giji (3) 

where gjj = gj#jj is the bare GF for H , % 6 {1 = 2,2 = 1}. To take into account the 
self-energy parts from the e-e correlations and e-p scattering if/, o"^-, we obtain the self- 
energies Jjij = S° + (Tij. If there is no coupling between ipi and ip 2 other than tunneling, 
then cr 12 = CTu and cr 2 i = cr 2 2- 

The Dyson equations @ can be solved easily: G^f = [(gj' a ) _1 — S[' a ] _1 , Gf { ~ = 
G^£+~G^, where we have used the total self-energy X, = + 5Xj with 



= er r t 



tjaCTj. = a x ti 2 git\ 2 a x . (4) 



Here we have used five different self-energies. Their physical meaning is following: as de- 
scribed above, £? is the self-energy due to tunneling to the neighboring lead; er^ is the 



self-energy due to interactions Hj without tunneling between quantum dots; S2i is the 
self-energy due to tunneling to another dot (Note this self-energy is calculated using the 
renormalized GF gj); finally, Sj is the total self-energy. 

From the Dyson equations, the current can easily be calculated at different dots: 
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where the superscripts (a,j) E {(L, 1), (R, 2)} denotes the current tunneling from the lead 
to the dot-i, and (j) is the inter-dot current calculated at dot-j. 

II. SIMPLIFIED CURRENT EQUATIONS: APPROXIMATIONS 



In general, the introduction of <Tij makes the current in Eq.(g) non-conserved ||12|| . To 
obtain a current obeying detailed balance, current- conserving approximations for <x are 
required. However, since we adopt HFA here, when Hj is neglected then the current in 



Eq.(pf) is automatically conserved [11]. The current is now written in a simpler form: 



J = - j duj(f L u - f R u )T r G^G^t^r^tL 



(6) 



where f[ R is the Fermi-Dirac distribution function and gj is defined in Eq.(^j). Recall that 
ti2gir 2 git{ 2 is just the imaginary part of self-energy S2i due to tunneling between quantum 
dots. The current equation @ is very similar to that for a single dot system, i.e. replace 
r a due to tunneling to the lead by = (S2j) r — (82j) a due to tunneling from dot-i to 
another dot-j. 

The electron-electron (e-e) and electron-phonon (e-p) interactions have two kinds of 
effects: renormalizing due to 3?cr r,a and changing the distribution of electrons due 
to <t k (and ^scr r ' a ). In the former effect, we simply renormalize E^, so the current 



equation is the same as Eq.(^). After we include the higher order (than HFA) energy 
shift, Eq.(|6|) contains rich correlation effects. One example is for the two-impurity model 
in the Hartree approximation studied by Niu et al 0. One might attempt to calculate the 
interacting NGFs using higher order (than Hartree approximation) truncation of the EOM 
and then calculate tunneling currents using Eq.(^|). This is questionable for double dot 
systems except when the truncated EOM is in the form of a Dyson equation. Otherwise, 
one needs to re-derive the current equation using the truncated EOM. For large quantum 
dots with more than 40 electrons, only the quasiparticle states near the Fermi surface are 
relevant to the resonant tunneling. The spectral functions of these quasiparticle states do 
not change significantly with adding/subtracting a couple of electrons in the dots except 
for a slight energy level shift due to the "charging effects". This has been confirmed by 
microscopic calculations by Wang, Zhang, and Bishop p0| , who also found that the many- 
electron energy levels are approximately equally spaced for quantum dots with more than 
30 electrons. Thus for quantum dots with a large occupation numbers, the energy shift 
due to e-e interactions can be treated phenomenologically (see Sec. [TV]) and the level spacing 
can be extracted from experiments. The scattering-in/out effect due to cr K is difficult to 
treat correctly in various approaches due to the detailed current balance problem mentioned 
above. However, qualitative effects can still be obtained after simplifying Eq. @ appropriate 
to different conditions, in a similar fashion to our discussions of nonequilibrium distributions 
below. 

Up to now, our equations are general for multi-level dots. To study a specific system, 
we need to know the matrix elements [t^]^ (hopping) and [IY| H .. Let us first examine the 
relevance of various matrix elements. The elements [1^]^. for li ^ lj are due to tunneling 
from different levels and lj to the same state in the lead. These are usually not small 
compared to the diagonal ones, and they describe the mixing of levels in the dots due to 
tunneling to the leads. If we assume the energy of the states in the tunneling window 
(i.e. /ii < Ei < yU 2 ) is much smaller than the barrier height between leads and dots, then 
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[TiJ^. = Tj. The matrix elements [t^]^ describe wavefunction overlaps of different levels 
between two quantum dots. The value of these matrix elements are different but roughly 
of the same magnitude if their energies are comparable. Similarly, if we assume the energy 
of the states in the tunneling window (i.e. fii < Ei < /i 2 ) is much smaller than the barrier 
height between two dots, then we can use the following assumption [ti 2 ]y. = t (thick-shell 
model f9|). Using the assumption [Tj], .,. = Tj, 

1 3 



■ r r,al 



1 1 j 1 



k 



» J 



• 1-(E°)'.°£ fc -fer 



(7) 

* ^ J- 



• r ' a - ^-sS 8 )" 1 ] 



where u£'° = [(j 

In the experiments r i; |t| <C Ai? (spacing between energy levels in the dots). Using 
Eq.(^), it is easy to show that near the resonance to ~ e; : for ij = i Q , [gi]^. ~ l/^i; for 
h ? l 0] [g[] lih ~ l/[(/, - Z )A£]; for ^ J, ^ / , [g] Wj ~ (u, - Q )/[(/ l - J )fe - ^o)(A£) 2 ]; 
for /j = l , lj 7^ Z , [giJ^j ~ !/[((?' — ^o)A-E]. The properties of the [Gy^ are similar. Thus, 
at Vi, \t\ <C AE, only the tunnelings between levels close to each other are important. 

For the resonant peak structure, we can simplify Eq.(^) near resonance, /II > ~ 

(2) 

e\ > iir, by taking into account only the resonance level contribution: 

J * / Mh ~ /fi) (.- f g'-^J?(r 2 /2 + r 12 /2F 

e /• . /l-/r r 2 r,|t| 2 

aw- 



h J^{u - 4 1} ) 2 + ( r i/2) 2 ( w _ e p) _ |f ' 2( — )2 + (fiL + Ei 



|f| 2 



' 2 ( u _ e W)2 +(ri/2) a-' T ^2 T 2 ( u _ e W)2 +(ri/2) a/ 

where we have used definitions ir 12 = (5S 2 ) r — (<5S 2 ) a and 5e^ 2 - ) = 9?(<5E 2 ) r . Near resonance, 

^ (!) ( 2 ) ^ T l*P X (2) MV-eif) 

a*l > i ~ i > r 12 = Fi (1) a ' ; *e< > = (1) v 
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III. NONEQUILIBRIUM DISTRIBUTIONS IN THE DOTS AND POPULATION 

INVERSION 



Before discussing the resonant current, let us first consider the nonequilibrium distribu- 
tions of electrons in the two dots. The distribution at each level U of dot i is given by 



with A(e}; J )r\ = T,\t\ 2 Eki m m M [§iLm ^ e energy level broadening due to tunneling be- 
tween dots. When the level < ^r {^r < according to our convention), f[ D {^) = 
fji D (4f) = 1) ^2(^2) = 1> ^ ms ^ eve ^ m ^ ne * s occu Pi e d; when the level ej^ > jj, L , 
fL D ( e< h) = fR D ( e( h) = 0' ^2(^2 ) = 0' tms ^ eve ^ ^ s em Pty- ^From Eq.(||), the energy levels 
in these two cases do not contribute to the resonant tunneling. A more interesting case is 
for levels lying between /i^j and /i^, where f[ D = 1 and f^ D = 0, electrons are in nonequi- 
librium states, and F 2 2 ~ r^+x^i strongly depends on the energy spectrum in the first dot. 
Using the properties of [g[} l d ., when a \eve\-i u of dot-2 is in resonance with a level in dot-1, 
A2 ~ |t| 2 /(4rf), the occupation number of level-i u is F 2 u ~ ifTrv+FF' wnen a level-i^ of dot-2 
is not in resonance with levels in dot-1, A2 ~ (\t\/AE) 2 <C 1, the occupation number of 
level-id is given by Ff d ~ \t/AE\ 2 T 1 /T 2 < ri/r 2 . 

There is a population inversion in the dot-2 (|t| ~ T) when e iu > e\ . the population 
inversion in dot-2 can be achieved due to the resonance of a higher level-i u (relative to an 
active level-z^) in dot-2 with a level in dot-1. The nonequilibrium electron distribution in 
dot-1 is similar. When the level e^" < [ir, FlAeP ) = 1, this level in the dot is occupied; when 
the level > /i^, F^ie^) = 0, this level is empty. When /i R < < /i^, F^ ~ r ?l r ■ 
When this level is not in resonance with a level in dot-2, Ai ~ (\t\/AE) 2 <C 1, ~ 1; 
otherwise, Ai ~ |t| 2 /(4r|), F t \ ~ 1/(1 + |t| 2 /(4rir 2 )). Thus, with \t\ < r i; Ffc ~ 1 when 




i=l,h=h 



i=2,li=l2 



HFA 



HFA 



(9) 
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This population inversion is similar to that in quantum cascade lasers for superlattices 
|T5| and other proposals for double quantum wells [p5|| . However, the population inversions 
proposed previously for quantum wells are achieved without elastic resonant tunneling. For 



example, in the double quantum well systems P5(| , the upper level is attached to the lead 
with high chemical potential fih and with E u < fi h , so the upper level is occupied; the lower 
level (either in the same well or in another well) is attached to the lead with lower chemical 
potential /i/ with E\ > fii, thus it is empty. Since transverse momentum is conserved and 
can be neglected the quantized energy levels in the wells are mismatched, there is no elastic 
resonant tunneling, and the system is in quasi equilibrium with two chemical potentials /x/j 
and The electrons can only tunnel via emitting photons. In the population inversion we 
proposed above, the energy levels are not in quasi equilibrium with the leads. The population 
inversion is achieved through resonant tunneling. Comparing to the proposal in Ref. p5[ for 
double quantum wells, we believe that our proposal is easier to realize practically. 

In the above analysis, we have neglected the inelastic scattering effects due to e-e and 
e-p couplings. In real systems, the e-e and e-p scatterings have important effects on the 
nonequilibrium distribution of electrons in the dots. At high temperatures, the nonequi- 
librium distribution will be thermalized due to inelastic relaxation. Here we concentrate 
on low temperatures, when the dominant contribution to cr comes from higher order e-e 
correlations which will affect the NGF's. Taking this fact into account, the F-function with 
a finite cr becomes: 

r 2 + 72 + A(e^)ri 

where i~fj = cr^ — cr - with the assumption [jj]^^ = 72- Physically, 72 is just the average 
level broadening due to e-e scattering, cr^f is responsible for changing of the distribution 
functions. Note that a^f ~ 72- Here A(e^ ) has been changed due to the introduction of a. 
So at resonance of level-^, ~ '^ra+^Kr^+^^+iffi » while F l ~ (M^/tTa + 72)- If 
(°~22)t~i u ~ (?22)t d ~i d , it is easy to show that Ff u > Ff d . In the extreme case that (022)^ ~ 72 
and (ct2 2 )^,1„ ~ 0, the population inversion survives when \t\ 2 > £ME|±T^i(p 1 _|_ ^ l ) 2 . For 
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the e-e scattering in the typical experiments 7» ^ 1^ f2~0fl, so the population inversion can be 



achieved at \t\ ^> T. Consequently, the e-e correlations reduce the magnitude of the resonant 
current but do not qualitatively change the resonance peak structure. 

IV. NONEQUILIBRIUM CURRENTS AND PEAK SPLITTINGS 

A. Nonequilibrium Currents 

In this section, we discuss the recent nonequilibrium experiments of Ref. @ using the 
general equations derived above. In these experiments, Tj, \t\ <C AE, so we can neglect 
the contributions of energy levels far from resonance (i.e. for experiment ||, we only need 
take into account 4 relevant levels). We use Eq.flf]) to calculate the currents, which neglects 
the high order scattering-in(out) processes. This is a reasonable approximation for low 
temperature experiments of quantum dots with electron numbers on the order of 90 [20], as 
we discussed above. 

The oscillations of the nonlinear resonant tunneling current vs. gate voltage in experi- 



ment 0] can be readily explained if we take into account the charging effect [13,22] for the 
renormalized HF energy levels and [see Fig.|T|. In Fig.^| we plot the conductance of 
the current tunneling through the DQD vs. the gate voltage V g \ on dot-1 using a constant 
bias voltage of 280/iV. (Only one group of the peaks is shown). The current is calculated 
using Eq.(||) for \t\ = 0.2fteV; Ti = T 2 = 4.0/xeV, lO.OfieV; the level spacing in dot-1 and 
dot-2 as A£ (1) = 150/xeV and A£ (2) = 230/zeV, respectively. 

Our assignment of the peaks in Fig.(||]) is the same as that in Ref. 0. We assume in 
Fig.(§) that the spectrum in dot-2 is fixed while varying V a \. Using the energy level spacings 
given in Ref. [0] , the spacing between the first and second peaks is ~ 80fieV, and the spacing 
between the second and third peaks is ~ !50fieV. However, in the experimental I-V curves, 
the first spacing is larger than the second one. One possible reason is that the level spacing 
is not regular, which is suggested in Ref. H. Note that the gate voltage scale for these 
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spacings are ~ 900fiV and 2500/iV^ respectively. In Ref. |2| they are converted to energy 
scales as 70fieV and 200fieV. Then the level spacing between level-2 and level-3 in dot-1 
is barely half of the average level spacing. Another possibility is that this discrepancy can 
be due to the charging effect. As discussed in Sec. [TIT], without resonance between the two 
dots, the "active" levels in dot-1 are filled and the "active" levels in dot-2 are empty; with 
resonance, the resonance level in dot-2 has population ~ |t| 2 / (4riT2) while the resonance 
level in dot-1 has population ~ 1 — |t| 2 /(4Fir 2 ) (\t\ <C Tj). So the levels in dot-1 are shifted 
down by 5t\ oc |t | 2 / (riT2) and the levels in dot-2 are shifted up by 5t2 oc |t | 2 / (riT2) - The 
net effect is that the current peak positions are pushed upwardby 5V oc |t | 2 / (riT2) - Recall 
that the peak height for the resonant current is J oc Itl 2 /^]^). In the experimental I-V 
curve, the magnitude of the first peak is much smaller than that of the second and third 
peaks, so the second and third peaks are shifted up much more than the first peak (due to 
large charging effects). Thus the first peak spacing can be larger than the second one even 
if we assume the active energy levels are close to their average values. 

Another consequences of this charging effect is that the Lorentzian line shape of the 
current will be changed. For the experiments in , the peak should be distorted slightly to 
the right side due to charging effects alone (e.g. the high gate voltage side has a longer tail). 
However, couplings to the other levels will induce asymmetric shoulders as observed in the 
experiments. Note that there is no 4-a peak [see Fig.0. The absence of this 4-a peak was 
argued in to be due to rapid relaxation of level-4 to level-3. As we discussed in Sec .|ITT| that 
if the energy of level-4 is below /i^, then level-3 is already filled without coupling between 
level-4 and level-3, so even at high temperature, there should be no substantial relaxation 
from level-4 to level-3. Since the experiments are conducted at very low temperature (34mK), 
the argument that when level-4 is lower than there is large relaxation from level-4 to level- 
3 is questionable. From our Fig.|l], we can see that, since all the levels below fiL in dot-1 are 
filled, they is a large charging energy between 4-a peak and this group of peaks, so the 4-a 
peak is "missing". The authors in [Q] found that the amplitude of the resonance peak will 
stay constant with increasing bias voltage, which was argued to support the large relaxation 
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of level-4 to level-3. However, as we discussed in Sec.|TT[ the population of the lowest 
active level stays constant with increasing bias voltage V even without inelastic relaxations, 
the amplitude of the resonance will stay constant with increasing of bias voltage, so this 
phenomenon cannot tell if there is large relaxation. Thus, all the observations in Ref. |2| 
can be well explained by our calculations without the assumption of large inelastic relaxation 
at low temperature ~ 34mK. 

From our calculations the half widths of these peaks are not ~ T/2, due to coupling 
between multi-levels in each dot, and the increase of T will eventually destroy the resonance 
as shown in Fig. (121). 



B. Peak Splittings 

In a symmetrical geometry, Ti = T 2 = T and = e[ = e, the resonant tunneling 
current Eq.(^) can be simplified as following: 

J * {jdu{f L - f«) [{u _ e y_ {lt f l ! l r y A)? + TW - (11) 

It is readily seen that there will be a splitting in the resonance peak 2(|t| 2 — T 2 /4) 1 / 2 when 
\t\ > T/2. The ^-dependence of the resonant peak splitting is shown in Fig.([|). The 
physical meaning of this splitting is clear in the strong interdot coupling limit \t\ ^> T: 

r 2 /4 r 2 /4 



J^lj du{f L - f R ) 



(12) 



(w-e- |t|) 2 + r 2 /4 (cj-e+ |t|) 2 + r 2 /4_ 

which is just resonant tunneling through a two level single dot system. This is not surprising, 
since in the strong interdot coupling limit, the two dots can be treated as a single system 
(for these two resonant levels-this is not true for other non-resonant levels), two degenerate 
levels will be split by 2\t\ due to interdot coupling. However, if the dot-lead coupling strength 
is stronger, then each dot is coupled to its own reservoir, and then weakly coupled to each 
other. If the two reservoirs are out of equilibrium, then these two dots cannot be treated as 
a single system. From Eq. (|TTD , there will be no peak splitting in this regime [23|. Notice 
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that, for \t\ ^ T/2, the peak splitting predicted here is different from that expected from a 
naive argument based on degeneracy lifting ||. When \t\ ^> T, the splitting ~ 2\t\, and our 
result becomes equivalent to that in CBT. 



V. SUMMARY 

In summary, we have developed a microscopic theory for the equilibrium and nonequi- 
librium transport through a double quantum dot (DQD). The analytical expressions for 
resonant tunneling current peaks were derived in both linear and nonlinear conductance 
regimes. We showed that the nonequilibrium electron distribution has important effects on 
the resonant tunneling through DQD systems in the nonlinear regime, and it is necessary 
to be take into account these nonequilibrium distributions to interpret nonlinear tunneling 
experiments. We also find that multi-level coupling affects the width of resonant tunneling 
current peaks. Using the exact results for a non-interacting system, we obtained a tunneling 



peak splitting for a symmetrical double dot 2y \t\ 2 — (T/2) 2 when \t\ > T/2. Most inter- 
estingly, we predict that the population of electrons in the active levels of one dot can be 
inverted in the nonequilibrium regime by changing the gate voltage to make the upper level 
in the 2nd dot match with an active level in the 1st dot. This mechanism of population 
inversion could be used in a double quantum well resonant tunneling diode. The population 
inversion proposed here is different from that of QCLs |L5 |, which is based on photon-assisted 



tunneling in a superlattice |24| or a double quantum well f25f . 
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FIGURES 
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(i) 

FIG. 1. Schematic energy potential landscape of the double quantum dots and leads. e\. s 
(i = 1,2) are the energy spectra in dot-1 and dot-2, respectively. The solid and empty dots 
represent populated and unpopulated states. a charging energy, produces populated states 

separated from the empty states. V g \ is a gate voltage applied to dot-1. V is a bias voltage. 



Fig.l of PRB, Zang et al 
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FIG. 2. Current through the DQD vs. gate voltage V g \ using a bias voltage V = 280/j.V. 
\t\ = 0.2fieV; = 150/xeF; AE ( ^ = 230/jeV 



Fig.2 of PRB, Zang et al 
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